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$^ , Abstract. In this short note we give incremental algorithms for the following lattice 

^l^i problems: finding a basis of a lattice, computing the successive minima, and determining 

^^ ' the orthogonal decomposition. We prove an upper bound for the number of update steps for 

every insertion order For the determination of the orthogonal decomposition we efficiently 

CO , implement an argument due to Kneser 



f^ ' 1- Introduction 

2 ' Many problems in computational geometry permit a natural computation by an incre- 

mental algorithm. Incremental algorithms process only one object at a time and insert it 
into a data structure. Most incremental algorithms follow an abstract framework: After 
processing a new object it is inserted into a data structure. It is first located where the data 
^ ' structure has to be changed {localization step). Then the data structure has to be updated 

^^ , locally {update step) in order to perform the insertion of a new object. 

Cn ■ Here we apply the incremental construction paradigm to the design of lattice algorithms. 

^V . Let Vi, . . . , v„i be vectors which span a Euclidean space E and let L ~ Zvi + • • ■ + Zvm 

^s ' be the lattice which is generated by these vectors. Suppose that we want to compute a 

\^ , property of L. First, we compute the property of the lattice Li = Zwi. Then we check 

^^ ' whether V2 G £i (localization step). If V2 £ Li, then nothing has to be done. If V2 ^ ii, 

fH , then we compute the property of the lattice L2 = Li + Zv2 (update step), etc. In every 

update step we compute a lattice basis for the new lattice Li which is computationally more 
expensive than the localization step. Hence, this algorithmic framework is attractive if the 
number of update steps is small. 
^ , After fixing notation in Section|2]we state an upper bound for the number of update steps 

for every insertion order in Section |3l In Section |3 we give algorithms for the following 
^ , lattice problems: an algorithm which finds a basis of a lattice given by a set of generators, 

C^ ' an algorithm for the computation of the successive minima of a lattice given by a complete 

set of generators (a generating set 5' of a lattice L is called complete if S contains every 
vector V ^ L \ {0} with ||w|| < raax^^s \\w\\), and an algorithm for determining the 
orthogonal decomposition of a lattice given by a complete set of generators. 

These considerations result in a simple meta algorithm with practical impact. It of- 
fers a significant performance benefit compared with straightforward implementations for 
classical algorithms. For experimental results see the technical report |,6|. This note is a 
concise version of this report where we in particular emphasize the incremental algorithmic 
framework. 
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2. Notation 

Let _E be a d-dimensional Euclidean space . Its inner product is denoted by ( • , • ) and the 
associated norm by || ■ || = y/{-, •). The d-dimensional unit ball is denoted by B^;. A point 
set L <Z E is called a lattice if there exist linearly independent vectors bi, . . . ,b„ G E 
such that L = Zbi + ■ • • + Z6„. Then, (61, ... , 6„) is called a basis of L and n is called 
the rank of L. The volume of i is given by volL — \ det(&i, . . . , 6„)|. A lattice L' C E 
is called a sublattice of L if L' C i. If the rank of L' and L is d, then the index of L' in 
i is [L : L'] = volL' / volL. The fc-f/z successive minima Xk{L) is the minimum value A 
such that XBd contains at least k linearly independent lattice points of L. We will need the 
following theorem of Minkowski (see e.g. Bit. 

Theorem 2.1. Let L C _E be a lattice of rank d. Then 

— volL< Xi{L)X2{L)---XaiL)yolBd < 2'^volL. 

3. Chains of Sublattices 

We want to construct a lattice L, which is generated by the vectors vi, . . . , Vm, incre- 
mentally. Update steps are necessary if Vi liVi + • • • + 'Evi-i, where i ^ I, . . . ,m. The 
next theorem gives an upper bound for the number of update steps. 

Theorem 3.1. Let vi, . . . ,Vm G ii^ be vectors which span E and which generate the 
lattice L. Define B — max \\vi\\. Consider the chain of lattices 

i — 1,...,771 



^m- 



(1) Zui C Zui+Zws C ... C Z-wi+ZwaH h Zw„ 

Then, in Q inequality holds at most d + log2{d\{B / Xi {L))'^) times. 

Proof. First we transform Q into a new chain of lattices which are all of full rank d. 
Choose indices 1 < ii < {2 < ■ ■ ■ < id < "m such that for all j S {1, . . . ,d} the 
rank of Zvi-^ + • • • + Zvi is j and ij is minimal with this property. We define the lattice 
L' — Zwjj + • • ■ + Zwj^ and consider the transformed chain 

(2) L' = TLvi + L' C Zwi + Zus + L' C . . . C Zt;i + Z-wz H h Zi;„j + L' = L. 

The number of inequalities in Q is at most d plus the number of inequalities in the 
chain (|2j. Define L'^ — Zui + • • ■ + "Lvi + L' , where i ~ 1, . . . , m. Since we have 

711 

yo\L'/yo\L=[L:L']=\{[L[:L[_^], 

2 = 2 

the number of inequalities in the chain (|2jl is at most the number of prime factors of 
vol L'/ veil/ which is at most log2(volL'/ voli). To finish the proof we apply Theo- 
rem l2.1l to the quotient volL'/ volL and use the fact Xd{L') < B. D 



An immediate consequence of Theorem B.ll is an upper bound for the size of a minimal 
generating set. 

Corollary 3.2. Let L C iJ be a lattice of full rank d. Let S* C L be a finite generating 
set of L. Then there exists a subset S' C S which generates L of size at most d + 
\og2{d\{B / Xi{L)Y) where B = max^,e5 ||i;||. 

Note that having long vectors in a lattice basis is not avoidable in general: Conway 
and Sloane |3| constructed a lattice in dimension 11 which is generated by its 24 shortest 
vectors but in which no set of 11 shortest vectors forms a basis. 
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4. Algorithms 

In this section we propose algorithms for lattice problems which take advantage of the 
incremental construction. The first two algorithms for computing a lattice basis and for 
computing the successive minima are straightforward. For the computation of the unique 
orthogonal decomposition we develop new ideas based on an argument of Kneser. 

4. 1 . Lattice Basis. For computing a lattice basis from a large set of generators we use 
an algorithm for computing a lattice basis from a small set of generators as a subroutine. 
Such an algorithm is the LLL algorithm for linearly dependent vectors of Pohst (see e.g. 
El, Chapter 2.6.4). Buchmann and Pohst (|1 1) showed for (a variant of) this algorithm 
that the number of needed arithmetic operations is 0{d + m)^ log B). For the incremental 
algorithm the number of arithmetic operations is linear in m. 



Algorithm 4.1 Lattice Basis 








Input: Generating system wi , . . 


• , V,n 


e E of the lattice L. 


Output: Basis 61, . . . , 6„ of L. 








n^O,L^{0}. 








for i = 1 to m do 








itvi^L then 








Use a subroutine to get n 


and a 


basis &i, . 


.,bnOfL + Zv^. 


L ^Zbi-\ h Zbn. 








end if 








end for 









4.2. Successive Minima. For computing the successive minima of a lattice our algorithm 
is similar to Algorithm 14. II However there are a few important differences: We need a 
complete generating system (see Section[0, the insertion order is no longer arbitrary, and 
in every update step it is enough to compute a basis of a subspace (instead of a lattice). 
Hence the number of update steps equals the rank of the lattice. 

Algorithm 4.2 Successive Minima 



Input: Complete generating system S — {v e L\{0} : 


||u|| < B}ofL. 


Output: Successive minima Ai(L), . . . , A„(L) of L. 




Choose V E S with minimal norm, S ^- S\{v}. 




n^ 1, L/^Kw, A„(L) ^ ||t;||. 




while 5 ^ do 




Choose u £ S* with minimal norm, S ^- 5'\{u}. 




if w ^ L then 




t/ ^ [/ + Rw. 




if dim U > n then 




n^n + 1, A„(i) ^ \\v\\. 




end if 




end if 




end while 





4.3. Orthogonal decomposition. A lattice is called decomposable if it can be written as 
an orthogonal direct sum of two non trivial sublattices. Eichler (|4|) proved that every 
lattice can be decomposed into indecomposable sublattices which are pairwise orthogonal 
and that the decomposition is unique up to order of summands. In [TJ Kneser gave a 
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constructive and much simpler proof. In this section we show how one can efficiently 
implement Kneser's argument. 

We are given a basis &i , . . . , 6„ G ii^ of the lattice L, a constant B, and a complete 
generating system S — {v ^ -^\{0} : \\v\\ < B}. We want to find the number of 
indecomposable sublattices r, indices ii — l<i2<---<ir<n<n + l = v+i and a 
basis b[, . . . , b',^ of L such that for every j £ {1, .... r} the vectors 6^ .,..., 6^ . _i form 
a basis of an indecomposable sublattice. 

Now we give Kneser's argument. 

Definition 4.3. A vector v G -^\{0} is called orthogonal decomposable if there exist 
x,y E L\{0} with v — x + y and {x, y) = 0. 

The orthogonal indecomposable vectors of S form the verticex set of an undirected 
graph G = {V,E). In G two vertices v^w E V are adjacent whenever (f,w) / 0. 
We decompose V into vertex sets Vi, . . . ,Vr of connected components of G. Then, the 
orthogonal decomposition of L is L = Li ± . . . _L L^ where Li is the lattice generated by 

Using standard algorithms from graph theory one can compute the connected compo- 
nents in time linear in 0(|y | + |i?|). In the following we show that in this case it is possible 
to compute the connected components in time linear in 0(|1/|). 

O'Meara observed in J8| that for the procedure above it is not necessary to determine all 
orthogonal indecomposable lattice vectors in S. The length decomposable lattice vectors 
are enough: 

Definition 4.4. A vector v G i\{0} is called length decomposable if there exist x,y E 
L\{0} with V = x + yand \\x\\ < \\v\\ and \\y\\ < \\v\\. 

On basis of this observation we propose the following algorithm. Its correctness follows 
from Proposition l4.6l In what follows we denote by vr^ the orthogonal projection of E onto 
the subspace spanned by Li. 

Algorithm 4.5 Orthogonal Decomposition of a Lattice 



Input: Complete generating system S = 


{v G A{0} 


■■h\ 


<B} 


ofL. 


Output: Indecomposable sublattices Li with L = Li 


_L .. 


LLr. 




Choose V e S with minimal norm, S ^- 


S\{v}. 








r ^ 1, Lr ^- TLv. 










while 5 7^ do 

Choose V E S with minimal norm, S 


- S\{v}. 








if^^El-ii^then 

Af ^ Zw + E,e 7 U. 
{ii,...,L,_|,,|}^{L, :*^ J}, 
end if 


-^■r-|,/| + l ^ 


- Af, 


r ^ r 


-|J| + 1. 


end while 











Proposition 4.6. At the end of each iteration the computed sublattices are indecomposable 
and pairwise orthogonal. 

Proof. By induction the sublattices Li, . . . ,Lr are indecomposable and pairwise orthog- 
onal. Let w be a shortest vector in S. If u ^ X]I=i ^»' '^hen v is not length decompos- 
able. In particular we have either Tri{v) — or 7ri(v) ^ Li where i = 1, . . . , r. Define 
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J = {j e {1, . . . ,r} : TTj{v) ^ 0}. One can choose vectors Vj G Lj, where j G J, 
which are not length decomposable and which are not orthogonal to v. In the graph G 
these vectors are all adjacent to v. Hence, Zu + J2ieJ ^J ^^ indecomposable and we get 

J2iei\J ^i -'- (^■^ + Y^jeJ ^j) because ■n^{v) = for i e I\J. U 

5. Acknowledgements 
We thank Martin Kneser for pointing out the reference to isl . 

References 

[1] J. Buchmann and M. Pohst, Computing a lattice basis from a system of generating vectors. Lecture Notes in 

Comput. Sci. 378, 54-63, Springer- Verlag, 1989. 
[2] H. Cohen, A course in computational algebraic number theory. Springer- Verlag, 1993. 
[3] J. H. Conway and N.J. A. Sloane, A lattice without a basis of minimal vectors. Mathematika 42 (1995), 175- 

177. 
[4] M. Eichler, Note zur Theorie der Kristallgitter, Math. Ann. 125 (1952), 51-55. 
[5] P.M. Gruber, C.G. Lekkerkerker, Geometry of numbers, North-Holland, 1987. 
[6] B. Hemkemeier and F. Vallentin, On the decompositon of lattices. Electronic Colloquium on Computation 

and Complexity TR98-52, 1998. 
[7] M. Kneser, Zur Theorie der Kristallgitter, Math. Ann. Ill (1954), 105-106. 
[8] O.T. O'Meara, On indecomposable quadratic forms, J. Reine Angew. Math. 317 (1980), 120-156. 

Universitat Dortmund, Fakultat fur Mathematik, 44221 Dortmund, Germany 
E-mail address: bhemkemeier@gmail . com 

Centrum voor Wiskundeen Informatica (CWI), Kruislaan413, 1098 SJ Amsterdam, The 
Netherlands 

E-mail address: frank . vallentin@gmail . com 



